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The instanton approach to spontaneous chiral symmetry breaking is reviewed, 
with emphasis on the connection to chiral random matrix theory. We extend the 
Ch ■ approach to discuss the finite density, zero-temperature behaviour of quark matter. 

^ ■ Since the instanton-induced interactions are attractive in both qq and qq channels, 

a competition ensues between phases of matter with condensation in either or both. 
qq ' It results in chiral symmetry restoration due to the onset of diquark condensation, 

a 'colour superconductor', at finite density. 

> . 

1 Why Use Instantons? 

The idea that the QCD partition function is dominated by instanton fluctu- 
ations of the gluon field, with quantum oscillations about them, has success- 
fully confronted many facts we know fifut the zero-temperature, zero-density 
• hadronic world (for reviews see refs. Era). Instantons have been reliably iden- 

tified in lattice simulations (for a review see ref. cl). From the theory side, the 
instanton vacuum constructed from the Feynman variational principleu, gives 
an example of how the necessary 'transmutation of dimensions' can actually 
happen in QCD, meaning that all dimensional quantities can be expressed 
y-H ' through the QCD scale parameter, Aqcd- 



2 Instanton Partition Function 



5h 

' The strategy of the instanton approach is to replace the partition function 

of full QCD with a sum over instantons and quantum fluctuations about 
them. In particular, the result is an integration over collective coordinates 
of instantons. A finite-density extension of this model is presumed valid for 
quark chemical potentials less than approximately 600 MeV, a scale set by 
the instantons themselves. 

The QCD partition function is approximated in two steps. First, the 
gluon degrees of freedom are replaced with a sum over an ensemble of in- 
stantons (Z's) and anti-instantons (J's) with quantum fluctuations above the 



trento: submitted to World Scientific on February 1, 2008 1 



non-perturbative vacuum. Explicitly, we write 

where are quantum fluctuations above the instanton vacuum and £ are 
the collective instanton coordinates comprised of positions, sizes, and colour 
orientations. With this ansatz, the partition function, including the quark 
chemical potential, is 



? = V 1 /jt)d -[rm.«(0 det(iff + im - ifi^)det{ip + iM) 

^N+\N-\J 5 p det(*y+iM)det(tp + im) ' 

Here m is a current quark mass, whereas M is a Pauli-Villars regulator mass. 
This expression provides a lower bound on the partition functiorO. 

Next, the quark determinant is partitioned into high and low momentum 
contributions. These regimes are divided by an arbitrary mass parameter 
Mi, which is chosen large enough such that Dethigh ifi factorizable and small 
enough so that Det/ oto is saturated by the zero modeaj. With it, we partition 
as 

Det = Det; OJi , • Dethigh 

dct(i'ft+ im — z/i7 4 )det(z^+ iMi) 



Dcti nn , — 



det(if+ iMi)det(ip+ im) 
det(if+ iM x )det(ip+ iM) 



Det '" 5 ' 1 det(ip+iMi)dct(if+iM) (3) 

If we consider the fields to be perturbative and subsequentially disregard 
their contributions, we may perform the averaging over the background gauge 
fields, represented by instantons. For the low momentum part one gets the 
forma: 

z low = y: Vtttt- J n **idpidOi d( PI ) n ( ^ ) < 4 . 



^ N+lN-l 

N± 



The T matrix is of size N+ x N-, where N± are the number of J's and J's, 
and is comprised of the overlaps of would-be zero modes, 

t ii = / d. 4 x <t (x - Zi,pi,Or,(j,)(ip- 2/174) $0(2; - z T ,pj, Oj;p) . (5) 



These integrals are not only functions of the instanton coordinates, but have 
parametric dependence on chemical potential and temperature as well. 
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3 Chiral Symmetry Breaking by Instantons 

The instanton partition function may be analyzed in three ways, each concep- 
tually quite different but leading to equivalent results. In each case, the main 
result is spontaneous breaking of chiral symmetry in the vacuum. This phe- 
nomena can be interpreted as a derealization of the "would-be" zero modes, 
induced by n the background instantons^resulting from quarks hopping be- 
tween themQ. It was first noticed in ref. □ that there is a far reaching analogy 
between chiral symmetry breaking in QCD and the problem of electrons in 
condensed matter systems with random impurities. The acquisition of a dy- 
namical (sometimes called constituent) mass by a quark is fully analogous to 
the appearance in the Green function of an electron in a metal of a finite 
relaxation time (but in our case this time depends on the momentum). The 
appearance of the massless pole in the pseudoscalar channel corresponding to 
the Goldstone pion is analogous to the formation of a diffusion mode in the 
density-density correlation function. For the recent development of these and 
related ideas, see Refs. and references therein. 



3.1 Random Matrices Approach 

Eq. (|J) is a particular (and historically first) example of a matrix approach 
to chiral symmetry breaking. The matrices here are made of the overlaps Tjj 
of the would-be fermion zero modes in the background of instantons whose 
positions, sizes and orientations are random. To make a quick estimate of the 
spectral density of the Dirac operator, v(X), one can use the cluster decom- 
position in the elements of the Tjj matrixH: 

u(X) = 5>(A-A„) = J ^ e"" A , (6) 



where the averaged trace over an ensemble of matrices is expanded as 



Tr e isT = iVexp <^ is— Tr T + — 
y \ N 2 



1 



— Tr T 2 T Tr T 

N N 2 



(7) 



The first term vanishes upon averaging. From the leading cummulant one has 

where N is the total number of I's and 7's and the parameter k is the average 
overlap of zero modes: 

= ^[TuTrA = const • yp 2 , (9) 
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Figure 1. Schematic eigenvalue distribution of the Dirac operator. The solid lines are the 
zero mode and free contributions, the dashed line an estimate of the full spectrum. 



with p 2 being the average size of instantons obtained from the size distribution 
d( Pl ) of eq.(|). 

Using the first cummulant in eq. (0) is, however, incomplete. Resumming 
all contributions yields the Wigner semicircle spectrum: 



"W = ^V 1 -^- (10) 

The chiral condensate is obtained directly from this result via the Banks- 
Casher relation. With the average instanton density N/V and sizes p com- 
monly used for the instanton medium, one gets the reasonable 

We note that the spectrum-smearing parameter k is a small number - approxi- 
mately 100 MeV. One likewise obtains a good approximation to the pion decay 
constant, findine_iL ~ 100 MeV, comparable to its well-known experimental 
value of 93 MeVOfl 

There is naturally more to the spectrum than the zero modes; mini- 
mally, the high-momentum modes which can be approximated by those of 
free quarks. The eigenvalue density obtained from this part of the spectrum 
is 

KA)«^A 3 . (12) 
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The full spectrum, shown schematically in Fig. hi is a combination of this, 
the Wigner semicircle, and the intermediate modes which are ignored in this 
treatment. _ 

It was noticed by SimonovEl that the partition function (^) can be gener- 
alized to the schematic form of a random matrix integrated with some weight 
whose precise form is not essential because of the 'universality'. This allows 
an economical rewriting: 

Zi ow = JdT P(TTt)det 

With an exponential weight, the form of eq. ( |13| ) is the partition function of 
chiral random .matrix theory, from which one obtains the Wigner semicircles, 
see also Ref. h. Later on, eq. ( Jl3| ) became a base of a now well-developed 
study of the microscopic spectral density of the Dirac operator and their 
correlations; see J. Verbaarschot, these proceedings. 

At finite density, the expressions must naturally he-modified. This has 
been done in a schematic way with random matriceal3. In the instanton 
picture, we retain the microscopic physics of the zero mode solution at finite 
chemical potential, which shall be discussed in what follows. 



im T 
T* im 



(13) 



3.2 Propagator Approach 

The random matrices approach is, unfortunately, limited to the spectral den- 
sity: it cannot be used to study the ^-dependent correlation functions related 
to chiral symmetry breaking. Knowing the microscopic origin of the random 
matrices Tjj as due to the overlap of the instanton zero modes, allows one to 
calculate those quantities. 

The simplest quantity is the ensemble-averaged quark propagator com- 
puted in refs. (second and third reference) and in ref. c3. The propagating 
quark interacts with the background via the would-be zero modes, quantified 
in the same overlap matrices T/j used previously. It has the averaged form: 



S(x -y) = §*{x - ZI ) . $j(y - zj) = . (14) 

im + lij p—iM(p) 

The result is a dynamical quark mass M(p) = const ■ y // p 2 N/VF(pp), where 
F{pp) is a function known analytically and features a natural cutoff at mo- 
mentum pp — 1 . 

Closing a quark loop with this effective propagator generates a 'gap equa- 
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tion' encoding the instanton density dependence of the dynamical massB, 

K-AN f d 'P CIV* 
V ~ C J (2tt)4 p2 + M(p )2 • ^> 

Since it has a solution with M(p) ^ 0, even for a zero current quark mass m, 
it implies spontaneous chiral symmetry breaking—. A nonzero current quark 
mass m can be also incorporated in this approachBllil. Knowing the propagator 
S(x, y\ m) one can find the determinant for low-momentum modes according 
to a general formula 

/•Mi 

Det; ou; = exp / dm' [S(x, x; <m') — Sq(x, x; m')] 

J m 

= exp|-A^y dm 'Y2 {^ m ' 2 + 4k2 ~ m ') | ■ ( 16 ) 

Comparing it with another representation, 

flf A 2 + m 2 \ 

Deti ow = exp I- J dXv(X) In y + M 2 J > (17) 

one recovers the semicircle spectrum, eq. jlQ). It demonstrates the mathe- 
matical equivalence of the propagator and the random matrices approaches, 
though the first method ris, more powerful as it allows one to compute also 
x-dependent correlationsErQ. In both cases the large-iV c approximation has 
been exploited. 



3.3 Effective Lagrangian Approach 

In the previous section the effective quark propagator was obtained by aver- 
aging the instanton background's effects on a single quark. One can recover 
identical physics by first averaging over instantons at the level of the partition 
function. The result is an effective action for quarks which contains instanton 
effects in induced multi-quark interactions. The would-be zero modes serve 
here as a bridge, passing information from the instanton vacuum to the ef- 
fective quarks through the induced vertex. The consequent interactions are 
vertices involving 2Nf quarks, commonly cited as 't Hooft interactions since 
he was the first to specify the proper quantum numbers. 

The overall strength of this interaction is not fixed once and forever but 
rather has to be integrated over all possible values, since it technically appears 
as a Lagrange multipliertl. Fortunately, in the thermodynamic limit one finds 
the integral is completely determined by its saddle-point result. Denoting the 
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co^phng constants A±, the result for the low-momentum partition function 
is 



Z = J d\+d\_ J Difj exp J J d 4 x ^(ip- ifrf^ip + \+Y+ f + \_Y Nj 



Y .(ln^t.-l) +JV _(ln^-L) } 



The left-handed vertex, 

=Jdu /n t 1 ^^^ 1 ] c^(E(p/ - n H ,u i7 

■4/a /i /P/)^/>M)i / / fe ^^„ /0 ^t(fc /) _ M )^^/(fc /) ] (19) 

contains the matrix form factors 

•^O,/") = {P + <f(P, V) + > (p, —fi) = <p*(p, — /i)~(p + • (20) 

The tp(p, /j,) are the Fourier-transformed zero mode solutions and we use the 
notation notation a; ± = x^a^, where the 2x2 matrices — (±ic?, 1) de- 
compose the Dirac matrices into chiral components, and it is understood that 
fi written as a four- vector is [i a — (0,/i). A similarly defined Y~ carries 
right-handed quarks. All calculations discussed here are in the topologically 
neutral case, where N + = A_ = N and hence A = A + = A_. 

With a Fierz decomposition, the instanton-induced interaction can easily 
be made to resemble the Nambu-Jona-Lasinio model. However, in the steps 
leading to this result, microscopic information has been retained. Instantons 
not only provide a natural cutoff through the form factors and a dynamically 
determined coupling strength, but also account for anomalously broken U^(l) 
symmetry and correctly reproduce SU(4) Pauli-Giirsey symmetry when iV c = 
2. It should be also added that the naive addition of a nonzero current quark 
mass to the NJL Lagrangian fails to reproduce several known low-energy 
Ward identities, as well as the phenomenologically- known coefficients in the 
Gasser-Leutwyler chiral Lagrangian (the terms containing m 2 and m-p 2 ). The 
microscopic instanton approach preserving all symmetries of QCD is capable 
to correctly incorporate nonzero quark masses, and it does so in a rather 
nontrivial waj£3. 
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4 Competition Between qq and qq Channels 

Since the instanton-induced interactions ( [l9| ) support both qq and qq conden- 
sation, it is necessary to consider the two competing channels simultaneously. 
This means that one must calculate both the normal (S) and anomalous (F) 
quark Green functions. A colour/flavour/spin ansatz compatible with the 
possibility of chiral and colour symmetry breaking is 

(V> /m '(p)^>)> = SffiSiip)} for a, (3 =1,2, 

ty /a< (p)< /3J (p)> = 5 f g 8 a p S 2 {p)) for a, > 2 , 

(^ ai (p)^ W (-p)> = (Tp f R ai ti>)< j (-p)) = ef'<*W&F(p) , (21) 

where [7] refers to some generalized direction(s) in colour space, and it is 
this set of N c — 2 indices which signals the breaking of colour symmetry. 
In the particular case of N c — 3, where the colour symmetry is broken as 
SU(3) — > SU(2) x £7(1) and our ansatz considers the 3 channel, we will by 
convention take [7] = 3; for N c — 4 one can take [7] = 34 and so forth. In the 
event of colour symmetry breaking, the standard propagators (and ensuing 
condensates) will lose their colour degeneracy and the separation of S(p) into 
S\(p) and S2(p) becomes necessary; otherwise the Schwinger-Dyson-Gorkov 
equations do not closetl 

Written in the chiral L,R basis, the 4x4 propagators Si^ip) are of the 
form: 



S(p) 



G( P )i z( P )s ( P y 

Z(p)S (p)- G(p)l 



(22) 



Here the off-diagonal, bare propagator So^) 1 * 1 = [{p + i s modified by 

the scalar functions Z 12 {p), and is augmented on the diagonal by the scalar 
Gipip) which if nonzero break chiral symmetry. 

Using the instanton-induced interaction ( |l9| ) directly, without Fierz rear- 
rangement, one can build a systematic expansion for the F, G Green functions 
in the parameters X/N c and p/R. In the leading order in both parameters we 
restrict ourselves to the one-loop approximation shown in Fig. |^. It cor- 
responds to a set of self-consistent Schwinger-Dyson-Gorkov equations. An 
important /i-dependence enters through the form factors in eq. (|l9|). 

4-1 The System of Gap Equations 

This system of equations reduces to a set of three gap equations which self- 
consistently specify three condensates: two chiral and one diquark. They are 
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Figure 2. Schwinger-Dyson-Gorkov diagrams to first order in A. 



defined as the closed loops: 

9i 2 = —5^ — [ t^t A(k, u)G x a(Jfe) ; / = — =^ — / S^j B(k, u)F(k) . 
y ' 2V|-1 7 (2tt) 4 V ' w MV J N*-1J (2tt) 4 V ' w V 7 

(23) 

The functions A(fc, jtt) and i?(fc, ^t) are scalar functioios_.which arise from spin- 
averaging the matrix form factors, as defined in ref. 113. The quantities 
and A are linear combinations of the condensates 51,2 and /: 

Mi = (5 - A/N c ) 9l + {2N C - 5 + 2/N c ) g 2 , 

M 2 = 2 (2 - l/N c ) 9l + 2(N C - 2)g 2 , A = (1 + 1/N C ) f . (24) 

The Mi ,2 are measures of chiral symmetry breaking, which act as an effective 
mass modifying the standard quark propagation. On the other hand the 
diquark loop 2A plays the role of twice the single-quark energy gap formed 
around the Fermi surface. 

The solution of the gap equations depends on the vertex coupling con- 
stant, A, which itself is determined by balancing the instanton background 
with the condensates through its saddle-point value. This minimization of 
the partition function (ft8|) leads to 



AT A(N 2 — 11 

- = \(Y+ + Y-) = -L-JL 1 [ 2 . 9l Mi + (N c - 2)g 2 M 2 + 4/ A] . (25) 

This joins the gap equations to close a system of equations, numerically solv- 
able. 

Once this is done, the chiral condensate proper may be computed as an 
integral over the resummed propagator: 

- (H)Mink = i{i>H) E uci = 4 [2Gi(p) + (N c - 2)G 2 ( P )] . (26) 
4-2 Thermodynamic Competition 

For any given chemical potential, multiple solutions can be obtained for the 
gaps. These correspond to different phases of quark matter, and they are 
summarized as follows: 
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Figure 3. Left panel: Condensates for N c = 3 as a function of fj,. Shown are the effective 
quark mass M, the quark condensate —(■tpip) 1 ^ 3 , and the diquark energy gap per quark A. 
Right Panel: The quark density n q vs. /i. 



(0) Free massless quarks: g\ = gi = / = 0. 

(1) Pure chiral symmetry breaking: g\ = gi ^ 0, / = 0. 

(2) Pure diquark condensation: gi = g2 = 0, / 7^ 0. 

(3) Mixed symmetry breaking: g\ 7^ 52 7^ 0, / 7^ 0. 

The free energy, calculated to first order in A, is minimized in order to 
resolve the stable solution. The phase corresponding to the lowest coupling A 
is the thermodynamically favouredE2l. 

No solutions were found matching Phase (0), and the Phase (3) solution 
obtained disappears at relatively low chemical potential (/z » 80 MeV) and is 
never thermodynamically competitiveca. The remaining phase competition is 
then between Phases (1) and (2). In the vacuum, where \i = 0, one finds Phase 
(1) preferred - this is the standard picture. However, at a critical chemical 
potential /u c , defined by the ratio of superconductive gap to chiral effective 



M(/x c ) y8(JV c -l) 4 ' [ 1 

a first-order phase transition occurs. With the standard instanton parameters 
N/V = 1 far 4 and p = 0.33 fm, we find fx c ~ 340 MeV. This and the first- 
order nature of the phase transition are clearly seen in Fig. ^[ 

The quark density is physically more relevant than the chemical potential. 
Thus we calculate 

d 4 xj 4 (x) = -i J -j^-^ n(p) , (28) 
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Figure 4. Occupation number n(p) vs. p for Phase (1) (left panel) and Phase (2) (right 
panel) for fi = l/p = 600 MeV. 



where we have defined the occupation number density for quarks as 
as an intermediate step. The correct, conserved form of the quark current is 
given in ref. E3. Here we present only numerical results, in Figs. || and ^. 
In Phase (1), there is clearly an effective mass brought about by spontaneous 
symmetry breaking, indicated by the reduced Fermi radius. We stress that, 
despite the complicated four-momentum dependence of the interaction, the 
resulting occupation number density appears as a perfect Fermi step function. 
Cooper pairing, however, smears the Fermi surface, and this is evinced in the 
second plot. The residual discontinuity at \p\ — // is the contribution from the 
free, colour-3 quarks which do not participate in the diquark. Those that do 
are smoothly distributed as a Bosonic condensate. 

Performing the integration over three-momenta, we arrive at the density 
as a function of chemical potential. In the right panel of Fig. || this is shown 
for the equilibria states, demonstrating the large discontinuity at the phase 
transition, where the horizontal line signifies the quark density of equilib- 
rium nuclear matter. The phase transition occurs at an extremely low quark 
density, a^pni^ptual conundrum which is seeing ongoing discussion in the 
literatureBB33'LD. 

5 Conclusions 

We have discussed three ways of obtaining the spontaneous breaking of chiral 
symmetry in the instanton approach. While the physics of each seem super- 
ficially distinct, it has been demonstrated that they are equivalent. One of 
these approaches was extended for and applied to finite density quark mat- 
ter, where a chiral symmetry restoring, first order phase transition reorders 
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the quarks into a colour superconducting medium. These results agree with 
related works, but unfortunately also share a mysteriously low transition den- 
sity, a mere fraction of that corresponding to equilibrium nuclear matter. 
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